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The ability to maintain linearly polarized output in single mode 
fibers is essential for utilization of polarization dependent receiver 
circuitry. Our measurements with long lengths of fiber (200 m) indicate 
that we can find input polarization angles which yield essentially 
linearly polarized output. However, we found that these polarization 
effects are greatly influenced by the presence of physical stress on the 
fiber such as stress due to bending, twisting, mounting, and other 
variations in ambient conditions. We conducted several experiments 
on short length fibers where special precautions were adopted to assure 
repeatability of the measurements. Our results indicate the existence 
of a general theoretical model that predicts the output polarization 
characteristics as a function of input polarization and fiber length. The 
model assumes the presence of two asynchronous, orthogonal modes, 
uniformly coupled over the entire fiber length. The model, however, 
cannot distinguish between uniformly coupled and uncoupled mode 
cases based on the output radiation measurements. 



I. INTRODUCTION 

The polarization characteristics of "single" mode optical fibers have 
been the subject of several previous publications. 1,2 An understanding 
of the polarization sensitivity of such fibers is important in assessing the 
applicability of polarization dependent optical circuitry. 3 An additional 
implication of polarization sensitivity is the introduction of delay dis- 
tortion in "single" mode fibers. 

In this paper, we present a theoretical model based on the propagation 
of two orthogonal modes; it is shown that based on the measurement of 
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the ellipticity of the output radiation alone, we cannot distinguish be- 
tween the following cases: 

(i) The existence of two orthogonal modes, uncoupled with different 
propagation constrants. 

(ii) The existence of two orthogonal modes, with identical propagation 
constants and uniformly coupled by some means of periodic perturba- 
tion 

(Hi) The most general case — two orthogonal modes, uniformly cou- 
pled, but with nonidentical propagation constants. 

The organization of the paper is as follows. In Section II, we summarize 
the theory, leaving the details to the Appendices. Section III describes 
the detailed experimental procedure utilized to measure the radiation 
ellipse as well as details of each of the measurements. In Section IV the 
experimental data is compared with simple theoretical results developed 
in Section II. 

It should be noted that the assumption of synchronous uncoupled 
modes will not verify our data. An additional finding of importance is 
that the most general case indicates that for a fiber of any given length, 
excitation conditions exist at the input that result in linearly polarized 
output. 



II. THEORETICAL MODEL 

For completeness, we state the obvious: for synchronous uncoupled 
modes, the output radiation would be linearly polarized independent 
of the orientation of linear input polarization and fiber length. Experi- 
mentally, this is not observed and, hence, this model is ruled inappro- 
priate. 

As detailed in Appendices B and C, we assume spatially orthogonal 
modes whose propagation constants differ by A/?; we further assume that 
the modes are uniformly coupled with a constant coupling k. The fiber 
output radiation ellipse, 4 as described in Appendix A, would possess the 
following characteristics. If a and b are the amplitudes of the semi-major 
and minor axis components of the radiation ellipse, then their ratio is 

_ , a 1 ± [1 - sin 2 20' sin 2 2a] 1 / 2 
n-±7= (1) 

b sin 20' sin 2a 

and the orientation \p of the major axis of the ellipse is 

\J/ = 1/2 tan" 1 (tan 20' cos 2a) - rj/2 (2) 

where 0' =0 + tj/2 
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Fig. 1— Photomicrograph of single mode fibers, core diameter = 4.7 n, outside diameter 
= 133^. 



6 = input orientation of polarization with respect to x axis 
rj =tan~ ' 2k/A0 
«=1/2VAj8 2 + 4k 2 2 
2 = fiber length 

When R = or », the output is linearly polarized. From eq. (1), we see 
that this happens at B' = ±rmr/2, independent of fiber length. Thus, this 
model predicts that even in the presence of either phase asynchronism 
or mode coupling (or both), there are specific orientations of input po- 
larization, 

6 = -t//2 ± m7r/2 (3) 

for which linear polarizatio n is observ ed for all lengths. The measured 
relative phase shift 2a = v / A/3 2 + 4k 2 2, can be thought of as due to an 
effective A/3 e that includes the effect of coupling such that 

Aft, = VA/3 2 + 4k 2 (4) 

In our measurements, the reference angle B' = was always selected at 
first by rotating the input polarization to that angle for which the output 
was linearly polarized; then based on the observations of the radiation 
ellipse at the output, we cannot distinguish between the following cases: 
(i) asynchronous, uncoupled modes, (») synchronous, coupled modes, 
and {Hi) asynchronous, coupled modes. These cases are summarized in 
Table I. 

III. EXPERIMENTAL PROCEDURE 

Figure 1 shows a photomicrograph of the fiber when illuminated with 
white light. The fibers used had a core of pure Si0 2 and a cladding of 
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Fig. 2 — Electron micrograph of etched single mode fiber. 
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Fig. 3 — Experimental set-up. 



B 2 03-6Si0 2 . The core and the outer diameters were 4.7 Mm and 133 Mm, 
respectively. Due to diffusion of B2O3 in the MCVD fabrication process, 
the fiber core was unintentionally graded in refractive index and the Arc 
was lower than expected for these compositions. 5 Figure 2 shows an 
electron micrograph of the etched fiber. 6 The effective core-cladding 
index difference was An « 2 X 10~ 3 . Fibers of lengths usually less than 
3 meters were used in our short fiber length measurements. 

The experimental arrangement is shown in Fig. 3. Light from a po- 
larized He-Ne laser is passed through a polarization rotator to facilitate 
rotation of the input polarization at the input of the fiber. The light was 
chopped to provide a reference signal at the receiver and then focused 
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on the fiber by means of a lens. A small section of the fiber near the input 
end was immersed in glycerol to remove any undesired cladding mode 
that might have been excited. The output light was collimated by a lens 
and passed through a dan-Thomson analyzer before being detected and 
displayed in a PAR receiver. A Babinet-Soleil compensator could be in- 
serted and removed, as needed, between the collimating lens and the 
analyzer. 

Experimentally, we found the polarization effects on long and short 
length fibers were greatly influenced by the presence of any stress on the 
fiber and other changes in ambient conditions. To ensure repeatability 
of the measurement, it was necessary to take many precautions. The 
fibers were held as straight as possible by gently taping them on to a flat 
surface, at regular but not necessarily identical intervals to minimize 
inducing any stress in the fiber. The input end of the fiber was held by 
using a vacuum chuck except in one experiment where it was held by a 
mechanical holder. In that case, it was clamped first using as little 
pressure as possible and then taped. In all our measurements, the input 
polarization angle, which resulted in linearly polarized output, served 
as the input reference angle, i.e., 0' = 0. The analyzer position oriented 
to measure the cross polarized component served as the reference axis 
for the orientation of the output ellipse. If we alter the angle 6', the output 
in general will become elliptically polarized. In order to obtain the phase 
difference 5 between components parallel and perpendicular to the 
reference angle 6' = 0, we have to orient the Babinet-Soleil compensator 
parallel or perpendicular to the original reference axis of the analyzer. 
With the plunger of the compensator adjusted to obtain linear polar- 
ization at the output, from the plunger position and the calibration of 
the compensator, the phase shift 2a can then be calculated. In each 
measurement, where input polarization is rotated with a given fixed 
length of the fiber, or where the length of the fiber is varied keeping the 
input excitation angle 6' fixed, the output analyzer was oriented parallel 
to the major and minor axis of the output ellipse to obtain the power ratio 
of these components. Orientation of the minor axis was also determined 
from the angular position of the analyzer. 

IV. EXPERIMENTAL RESULTS 

4. 1 Fixed fiber length, input polarization varied 

In this experiment, an input reference angle was found such that the 
output polarization was linearly polarized; we then measured the ellip- 
ticity of the output radiation as a function of input polarization rotation 
about this reference. Figures 4 and 5 compare the experimental data with 
the theoretical results obtained by computing 20 logiofl from eq. (1) and 
evaluating the orientation of eq. (2). The best fit was for 2a = 59°, which 
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Fig. 4 — Measured and theoretical power output ratio, R, of the major to minor axis 
components" of output radiation ellipse as a function of input polarization angle (6'). Fiber 
length L = 2.39 meters. 



was also verified from the compensator measurement. However, this 
experiment cannot give a specified value for the phase retardation per 
unit length, i.e., 5/L( = Afi e ), since we can not determine explicitly the 
integral number of tt/2 phase shifts included in the entire length of the 
fiber. 



4.2 Fixed Input polarization, variable length 

4.2. 1 On axis excitation (d'= 0) 

In this experiment, after finding the reference axis, orientation and 
power ratio R of the output ellipse was measured; the fiber was shortened 
at the input end repeatedly by a small amount (~ 5 mm) and the ex- 
periment was repeated by reorienting the input polarization to obtain 
a linearly polarized output. From Fig. 6, we infer, as the theory predicts, 
essentially a linearly polarized output independent of the fiber length. 
The experimental limit of 38 dB indicated in the figure is the limitation 
imposed by the degree of polarization of the laser source. Our mea- 
surements on long lengths (-200 m) indicate the cross polarized com- 
ponent was down 32 dB. However, the output polarization was subject 
to severe variations due to physical stress and other changes in ambient 
conditions. 

The data (not shown) indicate the input polarizer angle and the output 
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Fig. 5 — Measured and estimated relative orientation of the ellipse ^ as a function of 
input polarization angle 0' for L = 2.39 meters. 



analyzer angle remained essentially unchanged with deviations, within 
experimental errors, indicating negligible angular rotation of the fiber 
with length. 

4.2.2 Slightly off-axis excitation (0'« 5° ) 

This experiment was conducted by offsetting the input polarization 
by an angle 6' (—6°). The input fiber end was held mechanically; the 
output end was successively cut and the measurements were repeated. 
As expected the output was elliptically polarized; as a function of length, 
the polarization changes from elliptical to linear and then back to el- 
liptical. From eq. (3) for a small angular offset, i.e., for small 6' 

l«l~^V (5) 



and from (4) 



sin 2a 



$ = - tan" 1 [26' cos 2a] 
~ d' cos 2a 



(6) 
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Fig. 6 — Essentially linear polarization is observed at the output when input end is 
shortened with on axis (0' = 0) excitation. 



Figure 7 shows the excellent agreement between theory and experiment. 
In fitting the \p data, it was necessary to include a "twist" term of the 
order of 0.16° /mm. This implied a periodicity of the twist to be about 
2.25 m, which is very close to the length of the fiber and was traced later 
to the holding and mounting arrangement at the input end of the fiber. 
When the fiber was allowed to lie flat and straight, and the input end was 
held using a vacuum chuck, the twist disappeared. 

From the compensator measurements after each cut, the average A/3 e 
was evaluated to be «0.0581 radians per mm. As seen from Fig. 7, linear 
polarization occurs at about every 54 mm and from eq. (5), it is obvious 
that this occurs when 2« = 0, ir, . . . etc. Therefore 2a = A/3L = -k with the 
result A/3 e = 0.0582 radians per mm, in agreement with the compensator 
measurements. This leads to an estimated value of effective index dif- 
ference 

An e = — A = 4" = 5.86 X 10" 6 
2tt 2L 

between the modes. Thus, it seems reasonable to state that the fibers 
cannot be easily fabricated to such tolerances and, hence, the polarization 
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Fig. 7 — Analytical predictions that R varies as (& sin 2a) -1 , and \p varies as (P cos 2a) 
as a function of length for small 6', confirmed experimentally. Period between successive 
linear polarization is 54 mm yielding a A/3 = 0.0582 radians/mm. 



problems observed with these experimental fibers will qualitatively be 
observed for all nominally circular, single mode fibers. 

4.2.3 Excitation at 6' = x/4 

If the fiber was excited at 45° from the reference angle 6' = 0, two 
modes of equal amplitudes would be excited; at the fiber output they 
would be out of phase by 2a = A/3(X- Therefore in this case, the polar- 
ization will change, as a function of length, from circular to linear and 
then back to circular with elliptical polarization in between. For & - tr/A, 
from eq. (1) and (2), 

\R\=—~, (9) 



and 



tan a 



yp = ±tt/4 



(10) 



In eq. (4), when cos 2a goes through zero, i.e., at a = ir/4 and \R\ = 1 
(circular polarization), \p changes sign. Experimental results compare 
very well with the analysis as shown in Fig. 8. 

V. CONCLUSIONS 

The fibers used in these experiments were nominally circularly sym- 
metric; however, the data on the elliptically polarized output radiation 
indicates that very small deviations from circular symmetry possibly 
exist and that these irregularities cause serious polarization effects. It 
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Fig. 8 — Behavior of output polarization with decreasing input length. R varies as (tan 
«) _1 and \p = ±45°. The beat period is the same as before — 54 mm. 



appears unlikely that the fiber fabrication process can be improved to 
the extent that such effects can be completely eliminated. 

The results indicate that input excitation condition^do exist for which 
linearly polarized output will be observed for a fiber of any length. 
However, the output polarization is subject to severe variations due to 
physical stress and other changes in ambient conditions. Our experi- 
ments can not determine whether any mode coupling exists. However, 
' if any mode coupling exists, it must be uniform, at least over the few 
meter lengths of the fiber measured. 

Polarization behavior of single mode fibers may hopefully be im- 
proved, by making the cores purposely elliptical, provided that any pe- 
riodic perturbation such as "twisting" of the core is sufficiently small 
and does not result in strong mode coupling. Large ellipticity, however, 
will increase the modal dispersion. Further study is necessary to char- 
acterize the polarization behavior of elliptical core fibers and, if they hold 
polarization, to arrive at an optimum degree of ellipticity that would 
provide a balance between polarization control and dispersion effects. 
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APPENDIX A 

Orientation and Ratio of Major to Minor Axis of Vibrational Ellipse of 
Elliptic ally Polarized Waves 

Consider two time-varying (S.H.) orthogonal electric field compo- 
nents 

E x = a ie >M+*i) (11) 

E y = a 2 eJ(° t+6 d (12) 

such that a? + a\ = 1. The resulting electric vector, as seen from Fig. 9, 
traces, in general, an ellipse. The difference in the effective propagation 
constants between the two waves is 

A(i= L^ = L (13) 

where L is the fiber length. Whenever a i or a 2 is zero or when 8 = rrnr, 
where m is an integer, the output wave is linearly polarized. Furthermore, 
when a\ = ai, and 8 = (2m + l)ir/2, the output is circularly polarized. 
The ratio R of the semi-major to semi-minor axis of the ellipse can be 
defined as 4 

i? = ^=tan x . (14) 

o 

where x is an auxiliary angle (— 7r/4 < x — t/4). By choosing an angle <f>(0 
< <t> < 7i72) such that 

tan = — (15) 

the angle \f/ of the orientation of the resultant ellipse, with respect to a 
reference axis ox, and the parameters and x are related by 4 

sin 2x = sin 20 sin 8 (16) 

tan 2\P = tan 20 cos 8 (17) 

In addition, 

al + al^aZ + b 2 (18) 

By eliminating x from equations (14) and (16), and with the use of eq. 
(15), we find 

a.a£, i**i-*j M*gl (19) 

b 2a\a.2 sin 8 

and from eqs. (17) and (15), 

f = \ tan" 1 [ 2 ° ia2 2 cos 8 1 (20) 

646 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1978 



2a, 




E x =a 1 ei^ t+& \ 



2a, 



Fig. 9 — Output radiation ellipse. 



APPENDIX B 

Relations Between Input Polarization Angle, Fiber Parameters, and Output 

Ellipse 

We assume a cartesian co-ordinate system, as shown in Fig. 10, with 
ox and oy coincident with the axes of slightly elliptical core or with that 
of the birefringent axes of the core of a single mode fiber as the case may 
be. Let us also assume that a linearly polarized plane wave is incident 
upon the fiber at an angle with respect to ox; therefore the amplitudes 
of the excited modes with orthogonal polarization are cos and sin and 
represent the x and y component respectively. 

Assuming a uniformly distributed coupling characterized by a coupling 
constant k between the modes whose propagation constant differ by A/3, 
the complex field amplitudes are given by 7,8 

E x = cos cos a — j cos (t; + 0) sin a (21) 

E y = sin cos a + j sin (7/ + 6) sin a (22) 

where the parameter a, proportional to the length of the fiber, is 



<x = V(AP/2) 2 + k 2 z . (23) 

and the parameter 77 which relate the degree of coupling and the asyn- 
chronism between modes is given by 

tan 77 = 2k/A/3 (24) 



POLARIZATION EFFECTS IN FIBERS 647 



0'= ir/2, fl=ir/2 - j?/2 
LP) 




0'=A0',(1/A0'Sin2o) 
(EP) 



Fig. 10 — Input excitation conditions and the behavior of the resulting output polar- 
ization. 



For reasons indicated in the text, we exclude the case when both k and 
A/3 are equal to zero. For the case k = and A/3 finite, rj equals zero. If A/3 
= and k is finite, rj equals tt/2. When both k and A/3 are finite, < i\ < 
t/2. 

By rewriting eqs. (21) and (22) similar to eqs. (11) and (12), we see 
that 



a\ = cos 2 9 cos 2 a + cos 2 (tj + 0) sin 2 a 
a\ = sin 2 6 cos 2 a + sin 2 (77 + 0) sin 2 a 



and 



a 2 — 0,2 = cos 20 cos 2 a + cos 2(77 + 0) sin 2 a 
The phase constant 8 is determined by 

2 sin(2 + 7?) tan a 



tan 8 = 



sin 20 - sin 2(0 + rj) tan 2 a 



(25) 
(26) 

(27) 
(28) 



Substituting (25) through (28) in (19) and (20) we find, after considerable 
trigonometric manipulation, that the ratio of the major to minor axis 
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of the resultant output ellipse is 



m l±fl- sin 2 (20 + 77) sin 2 2a] 1 / 2 
sin (20 + 77) sin 2a 

and its orientation 

\P = - tan" 1 [tan (20 + rj) cos 2a] - J (30) 

By appropriate choice of the sign, we see that |i? | in eq. (29) is bounded 
by 1 and °° or and 1. Whenever \R\ goes to or <*, we have linear po- 
larization. For a given length of the fiber z, therefore a, this occurs 
when 

= -- ± m - , m = 0, 1, 2, . . . (31) 

As r/ varies from to x/2, the input angle at which the output is linearly 
polarized varies from to — 7r/4. In the most general case, at this specific 
angle, phase difference 5, as seen from equation (28), equals either zero 
or multiples of it and is independent of a indicating that the two modes 
are either in or out of phase with each other at the output for all lengths. 
This clearly shows for a given fiber with a finite k and A/3, there is always 
an orientation of input polarization that would result in a linear polar- 
ization at the output. 

Note that in the special case 77 = 0, 5 = 2a and therefore, 5 is dependent 
on length; however, at = ±mir/2, as seen from (25) and (26), one of the 
components goes to zero. Thus the condition for input excitation angle 
to achieve linear polarization at the output, for all values of a, is given 
by (31) and holds good for all values of 77 in the range < \y\ < ir/2. 

APPENDIX C 

General Theoretical Model Used for Observation of Output Elllptlclty 

In this Appendix, we consider the general case illustrated in the pre- 
vious Appendix and show that two orthogonal uniformly coupled 
asynchronous modes can be resolved into a new set of orthogonal modes, 
but completely uncoupled. This implies, even if any coupling exists, as 
long as it is uniformly distributed, our results can be explained by as- 
suming uncoupled, asynchronous modes. Obviously, the measured A/3 
will now include the effects of coupling. 

If no mode coupling exists between the modes, i.e., k = 0, then from 
(23) and (24) 

a = ^z (32) 

2 
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and 



V = tan" 1 (j£) = (33) 



Under these conditions, from (21) and (22), 

E x = cos 6e-J a (34) 

E y = sin Oe + J« (35) 

We have shown in the previous Appendix that in the most general case 
< 77/2 < ±7r/2, there exists an input polarization angle such that = 
—77/2, for which the output polarization remains linear for all lengths. 
As shown in Fig. 10, if we rotate the coordinate system by —77/2, then the 
angle 0' in the new coordinate system is 

d' = 6 + 77/2 (36) 

But substituting (36) in (21) and (22), and writing the amplitudes along 
Ox' and Oy' as 

E' x = E x cos 77/2 - E y sin 77/2 (37) 

and 

E' y = E x sin 77/2 + Ey cos 77/2 (38) 

Using eqs. (34) through (38) we can easily show that 

E' x = cosd'e-i" (39) 

Ey = sin d'e+J" (40) 

Eqs. (39) and (40) are identical to (34) and (35) with being replaced by 
0'. Thus the coupled orthogonal modes can be easily resolved in terms 
of uncoupled orthogonal modes. Obviously, although a remains the same 
as in (23), it can be defined as in (32) to include the effects of the nonzero 
coupling as 



By substituting (36) in (29) we find the ratio of major to minor axis of 
output ellipse is now given by 

1 ± (1 - sin 2 20' sin 2 2a) 1 / 2 
sin 20 sin 2a 

Therefore excitation at the input polarization reference angle (0' = 0) 
for which R -*■ or °°, the output remains linearly polarized for all 
lengths. Then if we vary 0' at the input, eq. (42) represents the ratio of 
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axes of output ellipse. The orientation of the output ellipse is now given 
by 

4, = - tan" 1 [tan 26' cos 2a] - tj/2 (43) 
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